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We use the Bohr-Sommerfeld quantization approach in the context of constituent quark models. 
This method provides, for the Cornell potential, analytical formulae for the energy spectra which 
closely approximate numerical exact calculations performed with the Schrodinger or the spinless 
Salpeter equations. The Bohr-Sommerfeld quantization procedure can also be used to calculate 
other observables such as r.m.s. radius or wave function at the origin. Asymptotic dependence 
of these observables on quantum numbers are also obtained in the case of potentials which be- 
have asymptotically as a power-law. We discuss the constraints imposed by these formulae on the 
dynamics of the quark-antiquark interaction. 

o 

Q I. INTRODUCTION 

(N \ 

It is well known that nonrelativistic or semirelativistic constituent quark models can describe with surprising 
accuracy a large part of the meson and baryon properties. The mass spectra of hadrons have been intensively studied 
O ! since the pioneering works of, e.g., Eichten et al. Q, Stanley and Robson 0] and Godfrey and Isgur Q (see for 
■ example pi 03, E3 ) • Other observables, such as decay widths, have also been well reproduced within the framework 
' of potential models 0. The success of these calculations shows that it is possible to understand the major part 
of the hadron observables using the simple picture for the quark-quark or quark-antiquark interaction provided by 
t-H ' constituent quark models. 

In this paper, we will focus on the meson spectroscopy. Many potentials have been proposed to describe the 
meson properties 0, ||| . Despite of this diversity, the central part of the interaction used in usual models presents 
similar characteristics: an attractive short range part with a confining long range interaction. The prototype for these 
potentials is the so-called Cornell potential Q which can be used to describe a large body of the meson masses except 
| the masses of the pseudoscalar states for which a spin dependent and a flavor mixing interaction is necessary. Indeed, 
. the role of this additional interaction is too important, in this sector, to be neglected even in a first approximation. 
Recent lattice calculations confirm that the Cornell potential fits in a good approximation the static quark-antiquark 
"*^h , interaction. 

Most of the models found in Refs. |3,|(| have a phenomenological nature: whereas the general behavior of potentials 
are often dictated by physical considerations, the values of their parameters are most of the time obtained by trial 
and error from a comparison with some experimental data. The values of various observables which can be generated 
from potential models will be strongly constrained by this general behavior of potentials and by the values of their 
parameters; however, despite of the ease with which the numerical values of these observables can be obtained, the 
exact nature of these constraints is often very difficult to infer from the numerical results. 

It is the purpose of this work to show that the use of a Bohr-Sommerfeld quantization (BSQ) approach can provide 
unexpected insight into potential models. We use this method to obtain some general characteristics of constituent 
quark models. In particular, we will show that it is possible to obtain analytical expressions for the spectra, root 
mean square radii, decay widths, electromagnetic mass splittings or electric polarizabilities, which closely approximate 
the numerically exact results obtained from full quantum Schrodinger or spinless Salpeter approaches, and which can 
thus be used, e.g., to derive accurate predictions for the dependence of the observables on the quantum numbers. We 
will study the constraints imposed by these formulae on the dynamics of the quark-antiquark interaction within a 
potential model. 



II. THE BSQ METHOD 

The basic quantities in the BSQ approach [8( are the action variables, 

J s =Ip s dq s , (1) 

where s labels the degrees of freedom of the system, and where q s and p s are the coordinates and conjugate momenta; 
the integral is performed over one cycle of the motion. The action variables are quantized according to the prescription 



J s = (n s + c s ) h, 



(2) 
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where h is Planck's constant, n s (> 0) is an integral quantum number and c s is some real constant, which according 
to Langer should be taken equal to 1/2 @. 

To draw some general conclusions about potential models applied to meson spectroscopy we use the BSQ formalism, 
in the next sections, with the Cornell potential. The properties of various observables are representative of those which 
can be obtained from the usual constituent quark models. In particular, conclusions about asymptotic behaviors (large 
value of quantum numbers) of these quantities will be quite reliable since most of the potential models use a linear 
confinement. Moreover, we will show that we can also obtain information about asymptotic behaviors for any power- 
law confinement. 




A. Nonrelativistic calculations 

The nonrelativistic Hamiltonian corresponding to the Cornell potential reads in natural units (h = c = 1): 

(3) 

Of course since the interaction is central, the orbital angular momentum p^ = L is a constant of the motion. The 
expression for the radial momentum p r derives from the conservation of the total energy E of the system, 

p r =±- y/-2\ia r 3 + 2[iE r 2 + 2^nr- L 2 =±- \/g(r). (4) 
r ' r 

The radial motion takes place between two turning points, r_ and r+, which are the two positive zeros of g{r)\ the 
three roots, (k = 0, 1, 2), of g(r) are 



, , -P f9 + 2Trk\ E 
r k = 2 X /—co S (-^- + (5) 



where 



P = -— (E z + 3ok , cos 9 = 



3a 2 v ' 2 V /'''' 

1 



where 



. 2E 4 + 9anE - 27— L 2 . (6) 
27a 3 V / 

As ri < < r2 < ro, the turning points are r_ = T2 and r + = ro- 

Quantization of trivially gives L = I + c$; on the other hand, quantization of J r leads to the equation 

2/i r+ (En + 2k) K( V ) + 2fiEr+ (r+ - n) E( V ) - 3{£ + c$f U{ir/2, 7, rj) - 

3 

-W (n + c r )r +y /r + - r! y/2/j.a = 0, (7) 

y> — r_ T— 

rj=— and 7 = 1 , (8) 

r+ - n r + 

K(x), E(x) and U(Tr/2 7 x,y) are the complete elliptic integrals of the first, second and third kind, respectively 0, p. 
904], and n is the radial quantum number. This appears to be a rather complicated equation since it cannot be solved 
explicitly for the energy. However, it leads to very accurate results if we choose the Langer prescription c r — = 1/2, 
as it can be seen in Table H] where a selected set of masses, computed from Eq. J7J), are compared with the results of 
an exact calculation || . The accuracy of the results presented in Table [I] is representative of that obtained for the 
states which are not displayed here. One sees that, even for small quantum numbers, the errors introduced by the 
BSQ approximation are remarkably small; indeed the errors in the spectrum reported in Table^do not exceed 0.1%. 

One of the most useful byproducts of the BSQ method is to provide simple asymptotic expressions for the total 
energy for large values of £ and n. Indeed, for large angular momenta (£ ^> n), the orbits become almost circular, and 
thus r_ w r + and 77 sa 0. In this limit, 9 — n and all the elliptic integrals are equal to 7r/2; keeping the leading terms 
in Eq. Q, we obtain: 

E~U-) 1/3 t 2/3 (*»»)• (9) 
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This is a well-known result: a linear potential does not lead to linear Regge trajectories in a nonrelativistic description; 
correct Regge trajectories can only be obtained with confining potentials rising as r 2 / 3 [Tl|. 

For large radial quantum numbers [n 3> (), the orbits have a large eccentricity, and thus n ~ r_ and r\ ~ 1. 
Although i] = 1 cannot be directly inserted into Eq. Q) because of the singularity of the elliptic integrals, setting 
n = r_ = into Eq. and evaluating an elementary integral leads to: 

E ~(?f (|;f " 2/a <"»'>■ < 10 » 

This is a new result: the radial quantum number does not appear explicitly in the classical Hamiltonian, and thus 
"naive" semiclassical methods fail to reproduce correctly the spectrum in this sector 0,^3- Eqs. © and (|1U[) show 
that the spectrum has the same asymptotic behavior in t and n. The square of the ratio of the "slopes" , i?, assumes 
the value: 



r =[y) ' M 

that is, it does not depend on any physical parameter. In fact, these properties remain valid for any power-law 
confining potential. Indeed for V{r) ~ ar a (a > 0), the large-£ behavior is found to be: 

/ 2 \ V(«+2) /2 x a/(o+2) 

E ~fe) + f2 ° ,, ° +2, ( ' >>,1) ' (12) 

confirming the result of Ref. |llj |: for large n, the turning points behave as: 



r_ ~ 0, (13) 
and Eq. leads to an elementary integral which gives: 

/ a 2 \ V(«+2) / ^ \ 2a/(Q+2) 



^V(2?FJ ( " >>fl ' (14) 

where B(x,y) denotes the beta function 0, p. 948]. We will discuss in Sec. Ill Cl the implication of these results for 
the meson spectroscopy. 

The BSQ method can be used to compute other observables as well. In quantum theory, the value of an observable 
A is obtained from the wave function of the system as an average value. In the BSQ context, we have in contrast to 
evaluate an average over time according to 



(A) = i/ A{t)dt, 115) 



r 

T , 

where T is the period of the radial motion. As an example, we calculate the mean square radii of the states bound 
by the Cornell potential, which were reported in Table [I] Using the nonrelativistic equation of motion f — p r //i, we 
can write 

2x 2 /" 

V" 

and 



(0 = ^/ -dr, (16) 

T J r _ Pr 



r+ i 

T = 2u — dr. (17) 

Pr 



One obtains: 



(r 2 ) 



15a 2 



9an + 8E 2 + I AEn - 



3a(€ + c ) 2 \ K(ji) 



[i ) r\K{j]) + (r + — r\)E{r]) 



(18) 
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Inspection of Table[I]shows again that the accuracy of this formula is excellent. It can be used to derive nice asymptotic 
formulae; indeed one obtains 



(r 2 ) 



I \ 2 / 3 
(If I. 



1 \ 2/3 2 



£ 4 / 3 (*»n), (19) 

(r») ~ ( — ) = (3tt 4 ) 1/3 n 4 / 3 (n » £). (20) 

The behavior of the r.m.s. radius is the same for large £ and n. Comparison with Eqs. (|9I10|) shows that the r.m.s. 
radius becomes proportional at large I or n to the total energy. This behavior is a property of a linear potential. The 
ratio of the "slopes" of the relations I|19I2U[) does not depend on any physical parameter as it was the case for the 
ratio R for the energy. These properties are still valid for any power-law confining potential. For V(r) ~ ar a , the 
large-€ and large-n asymptotic formulae are: 



(r 2 ) 



( 1 



\aa/j, 



2/(a+2) 



£ 4/(a+2) (^ >n ) 5 (21) 



/ 2\ ( 1 V /(a+2) B{3/a, 1/2) / a* \ 4/( " +2) 



V2a^y B(l/a,l/2) V B ( 1 / a I 3/2 

These expressions can also be used to study the electric polarizability of mesons since the latter is proportional to 
(r 2 ) 2 i. 

Other observables can also be evaluated like, e.g., the various decay widths of the system. In the nonrelativistic 
reduction, the main ingredient in the calculation of these quantities is the modulus of the £ — wave function at the 
origin, ^(O)! 2 , a quantity which is not available in our formalism. However, use of the following relation valid in 
the nonrelativistic case 

l*(0)| a = £^V (23) 



2tt \ dr 

makes possible an analytical evaluation of the decay widths within the BSQ approximation. For a Cornell potential, 
using the definition (|15fl and the equation of motion f = p r j '/u, the square of the modulus of the wave function at the 
origin is found to be 

|*(0)| 2 = ^(a + K<l/r 2 », (24) 

with 

1\ H(7r/2,7,T/) 



r+(ri K(f]) + (r + - ri)E(f])) ' 

where 7 and 77 are defined by Eq. JHJ). In Tabled we show a comparison between this formula and values obtained with 
an exact calculation using parameters of Ref. J5J . The accuracy is here also remarkable especially for light mesons, 
with an error smaller than 1% for the nn states. The error is about 5% for the bb ground state. In this sector the 
Coulomb part of the interaction plays a more active role and more important errors are introduced since the BSQ 
result is only exact for a pure linear potential. However, the error for the ratio |\E , i5(0)/^ / 25(0)| 2 is always smaller 
than 1% even for the heavy mesons. 

It is also possible to derive an asymptotic behavior formula for ^(O)! 2 for a potential which behaves as V(r) ~ ar a : 

with the total energy E given by Eq. (|14|) . The asymptotic value for the ratio of the wave function at the origin, for 
two states with radial quantum numbers equal to m and n, reads 



* m (0) 



*„(0) 



/m\2(a-l)/(a+2) 

(n) ' (2?) 



It depends only on the value of a and the radial quantum numbers considered. 

Note that the BSQ approach can also provide analytical formulae for electric mass splittings because the wave 
function at the origin and the mean value of 1/r, which can be calculated with Eq. (|15fl . are the main ingredients for 
the evaluation of this quantity. 
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B. Semirelativistic calculations 



Similar calculations can be performed within a relativistic kinematics. The semirelativistic Hamiltonian correspond- 
ing to the nonrelativistic Hamiltonian reads 



i 2 

2 i i 9 K 



H = 2^p 2 + -§ + m 2 --+ar. (28) 

The orbital angular momentum p^ — L is still a constant of the motion. The expression of the radial momentum p r 
derives from the conservation of the total energy M of the system, 

p r = ± — Ja 2 r 4 - 2a Mr 3 + (M 2 - 2an - 4m 2 ) r 2 + 2Mn r + n 2 - AL 2 = ± ^- Jh{r), (29) 
2r 2r 

The polynomial h(r) is here of the fourth order. One can verify that h(r) reduces to g(r) in the nonrelativistic limit 
(rn — * oo). The radial motion takes place between two turning points, r_ and r + , which are two zeros of h(r); the 
four roots of h(r) are: 



''3 



1 / m nr-\ M 1 / /— r-—\ M , s 



with, 



A± = -{U + 2P)±^L 



= ~^( M +4ok + 8to 2 ), Q = — 



, -S (6\ 2P T / 27 

E/ = 2\ — cos - -— , cos# = - 



3 V.3/ 3 ' 2 V -S 3 

S = - W, T = -^P 3 + ^PR - Q 2 

R = ^-j (M 4 + 8Af 2 (ok - 2to 2 ) + 16a 2 (k 2 - 4L 2 )) . (31) 

The two turning points are r_ = r2 and r + = r 3 . 

Quantization of leads to L = i + c^; on the other hand quantization of J r gives the equation 

a x K(rj) + a 2 U(n/2, 7, 77) + a 3 n ( tt/2, —7, J? J + a 4 E(t}) 



-2na(n + c r ) ^/r 4 — r_ \/r + — r[ = 0, (32) 



where 



and where 



r4 — n r_i_ — r_ r + — ?*_ 

77= , 7 = , 33) 

T4 — r_ r + — ri r + — ri 

"i = -(r+-rx)(r- -n)(M + a(r 4 -ri)), 

«2 = — 2(r_ — r\) (an + 2m 2 ) , 
a 3 = -2a 2 r + r 4 (r_ - 77), 

«4 = — (r 4 -r_)(r + -ri). (34) 

This equation present some common features with the non-relativistic formula 0: it involves complete elliptic 
integrals, it cannot be solved explicitly for the energy and it leads to very accurate results if we choose the Langer 
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prescription c r — = 1/2. A comparison of the results obtained from this equation and exact calculations [j| is given 
in Table ITTI The accuracy is less good for the ground states in the relativistic case. The poorest result is obtained for 
the p state with an error of about 2.8%. But the convergence is quite rapid since the error already reduces to about 
0.2% for the p(1450) meson. 

We can also obtain some simple asymptotic expressions for the masses of the states for large values of I and n. The 
condition for circular orbits, r_ — r+, leads to 

2VU = \/A+ + x/Al. (35) 

Since the values of the total energy M and the angular momentum I are important, we have 

A+ = A_ =-{U + 2P). (36) 

These two last equations give U = —P from which we obtain an expression for cos(0/3). Comparison of this expression 
with the definition of cos 6 (|31|l imposes R = 0, and we find 

M-2V2\/^7 (^>n). (37) 

This is the expected result [il Il4l Il5j : a linear potential leads to linear Regge trajectories in the relativistic case. 

For large values of the radial quantum number, the large eccentricity of the orbits implies r*i ~ r_ ~ and 
r+ ~ J~4 ~ M/a. Setting n = r_ = into Eq. Q and evaluating an elementary integral leads to: 

M~2-vWon (n>£). (38) 

This result shows that the spectrum has the same asymptotic behavior in I and n. The square of the ratio of the 
slopes R assumes the value 

R = tt/2, (39) 

that is, this quantity is still parameter-independent. In fact, as in the nonrelativistic case, these properties remain 
valid for any power-law confining potential. For V(r) ~ ar a {a > 0), the large-^ behavior is found to be 



M^a l '^ +1 \a + l) {^j 
for large n, the turning points have the forms 



a/(a+l) 

{f>n); (40) 



r_ - 0, (41) 
and Eq. (JJJ leads to an elementary integral which gives: 



(42) 



To calculate the r.m.s. radius in the semirelativistic formulation we use the definition (|15f) with the equation of 
motion r = 4p r /(M — V(r)). The relativistic version of Eq. H16[) is found to be 

<r 2 ) = ^ [ + -(M-V(r))dr, (43) 

2T J r _ p r 



with 

T=ir (M -^^, (44) 

2 J r p r 

Using the expression (|29|l of the radial momentum p r the integration of Eqs. (|43I44|I leads to 

[r ' 3a2 faKiri) + p 5 E{rj) ' 1 ° J 
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with 

#l = -2Mn (M 2 + aK-4m 2 ) - 2M 2 n- 2a (n 2 -4(£ + c ) 2 ) 

+ a(M 2 + a« + 8m 2 ) (n(r + +n) -r_(r+ - n)), 

/3 2 = 24Mm 2 (r_ - n), 

/3 3 = -a(Af 2 +aK + 8m 2 )(r 4 -r_)(r+-ri), 

/? 4 = -2Afn -2k + o(n(r+ + n) -r_(r+ -n)), 

/3 5 = _ a ( r4 _ r _)( r+ _ ri ). (46) 

This rather complicated equation gives very accurate results as it can be seen in Table ITT1 where they are compared 
with exact calculations. It can be used to obtain some simple asymptotic formulae: 

(r 2 )~-l (£>n), (47) 



(r 2 )~—n (n»£). (48) 
da 

Like in the nonrelativistic calculations the asymptotic behavior is the same in i and n. Comparison with Eqs. 1371 - 
I38fl shows that the r.m.s. radius becomes proportional to the total energy for large values of quantum numbers, 
which is a characteristic of a linear potential. The ratio of the slopes which appear in these asymptotic formulae is 
parameter-independent. These properties are, here also, still valid for a power-law potential V(r) ~ ar a . In this case 
the asymptotic formulae are 

/ 9 \ 2/(«+l) 

(r 2 )~( — ) l 2 '^ (£»n), (49) 



1 /27r(a + l) 



2/(a+l) 



(r)~o — Ln ) («>£). (50) 



C. Discussion 



Application of the results obtained in Sec. Ill Al and III Bl to meson spectroscopy could prove to be of a great interest. 
Indeed, if we want to produce linear Regge trajectories (which are well observed experimentally in the light meson 
sector) within a potential model, we must use a potential which behaves at large r like r 2 / 3 in the nonrelativistic 
case, and like r for a semirelativistic kinematics. In both cases, we have just shown that the trajectories for radial 
excitation are then necessarily also linear; moreover, the square of the ratio of the slopes R is completely determined 
by the asymptotic behavior of the potential and the kinematics used, and takes here the value 

R=V3 (51) 

in the nonrelativistic case and 

R = -k/2 (52) 

in the semirelativistic case. This is an additional strong constraint, especially as these ratios are parameter- 
independent. But, since we have chosen the confining potential to reproduce the energy orbital trajectories (Regge 
trajectories), the asymptotic behavior of orbital and radial trajectories of other observables is determined. These 
results imply that the currently investigated potential models which use power-law confining interactions, can 

only describe a restricted class of experimental data. This remark could prove decisive in the (near?) future, when new 
experimental informations on the radial excitations of light mesons, which are still very scarce, become available. For 
example, if the experimental energy radial trajectories differ from a straight line, or if R does not assume the values 
of Eqs. H51|) or (|52l) . the understanding of the physics underlying the confinement could become more problematic, 
and a simple power-law potential would not be sufficient to describe the confinement of quarks (the above arguments 
remain valid if one considers the possibility that the asymptotic behavior of the Regge trajectory could not be exactly 
linear) . 
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If the experimental radial trajectories prove in the end to be linear, but if the ratio R differs significantly from the 
values of Eqs. (|5T)l or (JS^J, the introduction of a scalar component in the confining potential could be a possible way 
out; indeed we have performed calculations using the BSQ method, which show that this additional flexibility (which 
only makes sense within a relativistic approach) allows the ratio R to take any value between 7r/2 and 2. Of course, 
a quantum number-dependent confining potential could easily lead to a complete decoupling of the Regge and radial 
trajectories. A better kinematical treatment of the problem, as the use of full covariant equations (see for example 
Ref. ^(|)' could also alter the predicted large-€ and large-n behaviors of the trajectories. 

A limited set of additional experimental informations on the radial excitations of light mesons could already be 
sufficient to draw important conclusions. Indeed, the linear behavior of the Regge trajectories, which is governed by 
the long range part of the interaction, is already reached experimentally for the very first orbital excitations. The 
situation seems to be even more favorable for the radial trajectories, since, as the masses grow faster as a function of 
n than as a function of £, the asymptotic regime is expected to be attained at very low n; most potential models 0,EJ 
predict that this regime is already effective from the first radial excitation. 

To conclude this discussion, it is worth noting that the calculations presented in this paper allow to understand 
some well-observed properties of potential models. It is known that the use of a semirelativistic kinematics yields 
a better description of radial excitations such as p(1450) for the mesons and A^(1440) for the baryons. Actually a 
non-relativistic description gives in general too high masses for these states (see for example ^3 [SOU)- Indeed, 
correct description of Regge trajectories leads to heavier masses for radial excitations in nonrelativistic calculations 
since Eqs. <|51[1 and (|52[l show that the ratio R is higher when the Schrodinger equation is used. 

It is also known that experimentally the mass of the if* (892) meson is in good approximation the average of the 
masses of the p(770) and the 0(1020) mesons. This property is also verify for each orbital excitation of these states. 
This behavior is well reproduced within the usual potential models. But the strong relation between orbital and 
radial trajectories obtained for the energy in the previous Sections shows that this remarkable property will be also 
verify for radial excitations within a potential model description. This leads to a mass for the K*(2S) of about 1565 
MeV which is approximately the value found with the usual models; this value is just between the masses of the two 
possible candidates for this states if* (1410) and if* (1680). This is a major problem of usual potential models: the 
radial excitations of light strange mesons cannot be satisfactorily be described (see for example lalHI 20] ). 



In this section we show that a 1-dimensional BSQ approach can be used to derive some simple formulae which can 
be applied to the 3-dimensional I = states. For a central potential, odd states of the 1-dimensional Schrodinger 
equation remain solutions of the 3-dimensional i = Schrodinger equation (considering only the x > part of the 
x-axis). This property can be used within a BSQ approach. 

In a 3-dimensional calculation, if we use the Langer prescription, the centrifugal term never vanishes. But in a pure 
1-dimensional case, the absence of this centrifugal term simplifies the evaluation of the action variable and leads to 
very simple formulae. Indeed, in the non-relativistic case, we have 



With V(x) = a\x\ a , the quantization of the action variable J leads to an equation which can be solved for the energy: 



where we have changed n + 1/2 into (2n + 1) + 1/2 to only take into account odd states. This last result is an extension 
to the three-dimensional case of a formula obtained previously from a one- dimensional calculation |2lj | . This formula 
is very close from the Eq. I|14l) . but here it can be used, in principle, to approximate all the 3-dimcnsional I = states 
of a pure power-law potential. 

It is well known, from scaling properties (see for example J4|]), that the energy obtained with the Schrodinger 
equation for a power-law potential can be written as 



III. THE ONE-DIMENSIONAL BSQ APPROACH 



p = ±y/2fl{E-V(x)). 



(53) 




(54) 




(55) 



where e is a solution of a dimcnsionlcss Schrodinger equation. Since the BSQ approach gives the correct dimensional 
factor, the error depends only on a and n. Fig. n snows the evolution of errors with a for n = 1 and n — 2. For 
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the ground states the errors are about 1% (or smaller) for a < 4. As expected, the errors decrease rapidly with the 
increase of the radial quantum number. 

In general it is not obvious to understand why the formula H54[l works so well since the exact solution is not known. 
For the harmonic oscillator case this formula gives the correct position of the energy levels. But it is more instructive 
to consider the case of a linear potential. Indeed, for a = 1 Eq. I|54[) leads to 

f (- + 3/4)) 2/3 ^z 2 / 3 . (56) 

This expression is just the leading term of an expansion (z 3> 1) for the values of the zeros of the Airy function |22l 
p. 450]. The exact solution is obtained by the summation of all the terms present in the expansion. In this particular 
case we can see how the Eq. ()54|l approximates the exact solution and how evolves this approximation. The formula 
(|56[) gives much better results than variational methods using trial wave functions for which the accuracy decrease 
rapidly as a function of n (see for example [23, p. 267]). 

It is worth noting that a 3-dimcnsional calculation using BSQ approach introduces smaller errors than the 1- 
dimensional formulae obtained here. For example, for a linear potential the error is about 0.5% for the ground state 
when a 3-dimensional calculation is performed, and about 0.75% with the 1-dimensional formula. But in general the 
3-dimensional calculations do not lead to analytical formulae or lead to complicated ones. 

We can also calculate the mean square radius for a power-law potential. Using the definition l|15|) we can write: 



m-UsJ Mm\m^m> {n+m ) ■ (57) 

The calculation of the values of the wave functions at the origin leads to 

but this time with the expression l|54() for the energy. The evolution of the accuracy of formulae l|57|) and (|58|l with 
a and n is similar to that shown in Fig. ^ for the energy. 

The formula for the ratio of the wave function at the origin is given by 



,T, m\ 2 / i o U \ 2(a-l)/(a+2) 

y > m (0) f m + 3/4 x 



*„(0) 



(59) 



+ 3/4 

Some connection with previous general results can be done. If we suppose that m = 2 and n = 1, formula l|59|l shows 
that the ratio is greater than 1 if a > 1 (d 2 V(r)/dr 2 > 0) and smaller than 1 if a < 1 (d 2 V(r)/dr 2 < 0). This behavior 
is predicted by a general result obtained in Ref. [24| . Eq. I|54|l shows that the energy behaves with the reduced mass 
as E oc n~ a /( a+2 \ Thus we can calculate that 

^fW) K ^^/M, (60) 
d/i \ fi J a + 2 

This quantity is positive if a < 1 (d 2 V(r)/dr 2 < and dV(r)/dr > 0) and negative if a > 1. This property is proved 
for n = 1 in Ref. H|. 

The same calculations, for the energy and the r.m.s. radius, can be done for a semirelativistic kinematics but unfor- 
tunately they lead to complicated equations involving hypergeometric functions. The relativistic version of Eq. i|54|) 
can only be solved explicitly for the energy when one considers asymptotic behaviors. Thus further information about 
the observables cannot be extracted within the frame of a 1-dimensional semirelativistic approach. 

To conclude this Section, we show that the 1-dimensional BSQ approach can also be used to derive a simple formula 
giving the number of I — bound states (a generalization to I ^ states can easily be done). We consider an attractive 
potential which vanishes at infinity (for example a Gaussian or a Yukawa potential). We calculate the value of the 
radial quantum number n for E = 0: 

n= 2 ^ T ^W)dr-\. (61) 
t Jo 2 

In general this number is not an integer except if the energy level E = is a real solution. But the integer part of 
n gives the radial quantum number n max of the highest energy level. The formula for the number of bound states, 
N = n max + 1, of a given potential reads 



N 



(62) 
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where [x] denotes the integer part of x. For example, if V(r) is of the form V(r) = —af(br), we have 



N 



irb 



(63) 



The remaining integral is a pure number and we can see immediately the dependence of the number of bound states 
on the potential parameters. 

The same calculation can be performed in the semirelativistic case and we find 



N = 



y/V(r) 2 - 4mV(r) dr + 



(64) 



In this case the simple reduction performed in Eq. I|63[) cannot be obtained. 



IV. SUMMARY 



We have shown in this paper that the Bohr-Sommerfeld quantization procedure is an accurate and powerful method. 
It makes possible the derivation of reliable analytical formulae, from which one can easily study, e.g., the dependence 
of some physical quantities on the quantum numbers. Our study emphasizes the strong connection existing between 
the Regge and radial trajectories for the energy within a potential model, which could play an important role in 
elucidating the confining properties of the quark-antiquark interaction when additional experimental information is 
obtained on the radial excitations of light mesons; the study of other observables (decay widths, mean square radii, 
electromagnetic mass splittings,...) could put supplementary constraints on this interaction. 

We have also shown that even if semirelativistic calculations yield results quantitatively different compared with 
nonrelativistic calculations some common general features of observables are observed. The asymptotic behavior of 
orbital and radial trajectories, for the energy and the r.m.s. radius, are the same for a potential which behaves 
asymptotically as a power-law. Moreover the ratio of the slopes of these trajectories depends only on the value of the 
power of the confining potential. 

At last, we have shown in Sec. Ill Cl that 1-dimensional calculations lead to very simple approximated formula for 
the energy, the r.m.s. radius and the wave function at the origin in the £ = sector. We have also obtained formulae 
which give the number of bound states of a given potential for the nonrelativistic and semirelativistic case. 

Of course, a BSQ approach cannot replace a correct quantum description since it is only an approximate method 
which is not completely self-consistent. For example, we need to use the Schrodinger equation to give a definition of 
the wave function at the origin to be able to calculate this quantity with this formalism. Moreover some quantum 
concepts have no meaning in this framework. But we emphasize that for some cases the evaluation of observables 
(evaluation of average values) can be performed using this old method; even if calculations yield complicated formulae 
it is often possible to extract interesting information from an analytical relation and in this way guide full quantum 
calculations. 
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TABLE I: Energies, r.m.s. radii and wave functions at the origin for a selected set of mesons, obtained from the Schrodinger 
equation and from the BSQ approximation for the Cornell potential, using the parameters of Ref. §§ (masses are given in MeV, 
r.m.s. radii in GeV -1 , and wave functions at the origin in GeV 3 ). 
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TABLE II: Energies and r.m.s. radii for a selected set of mesons, obtained from the spinless Salpeter equation and from the 
BSQ approximation for the Cornell potential, using the parameters of Ref. |fj (masses are given in MeV and r.m.s. radii in 
GeV" 1 ). 
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FIG. 1: Evolution of the errors introduced by Eq. 1541 as a function of a for the ground state and the first radial excitation. 



